Topological charge pumping effect by the magnetization dynamics on the Surface of 
Three-Dimensional Topological Insulators 
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We discuss a current dynamics on the surface of a 3-dimensional topological insulator induced by magne- 
tization precession of an attached ferromagnet. It is found that the magnetization dynamics generates a direct 
charge current when the precession axis is within the surface plane. This rectification effect is due to a quantum 
anomaly and is topologically protected. The robustness of the rectification effect against first-varying exchange 
field and impurities is confirmed by explicit calculation. 
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I. INTRODUCTION 

A topologically-classified phase has universal properties 
which are robust against perturbation. Topological aspect 
emerges regardless of classical and quantum systems, and it 
is expected to be a useful tool to develop devices using non- 
trivial quantum phenomena. Recently, a topological insula- 
tor (TI), which is a gapped insulator in bulk and has a topo- 
logically characterized gapless edge state, has been drawing 
attention 1 - and the various exotic phenomena have been re- 
ported, such as the surface quantum Hall effect and the real- 
ization of the monopole dynamics^ For a TI in three space 
dimensions (3D), Qi et. al. argued that the surface state is 
described by the massive (gapped) (2+l)-dimensional Dirac 
fermion and non-trivial topological magnetoelectric effect ap- 
pears when the exchange field is present.— The low energy 
response to the electromagnetic field was shown to be un- 
derstood by the topological mass (TM) term (Chern-Simons 
term) 7 .—: 
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where d?x = dtd 2 x, m is the mass of the Dirac fermion and 
A is the £7(1) gauge field. 

One of the interesting aspects of TM is the quantum 
anomaly, where quantum fluctuations break the symmetry in 
the original Lagrangian. In fact, TM breaks the parity sym- 
metry which exists when Dirac mass is zero: the so-called 
parity anomaly. Originally, the quantum anomaly was dis- 
cussed in the context of the chiral anomaly in the decay of the 
7T meson, which has been experimentally well confirmedji 2 *^ 
Generally, the anomaly emerges from the regularization of the 
ultraviolet divergent diagram in the diagrammatic expansion. 
In condensed matter physics, it is known that the graphene 
exhibits exotic physics because of TM, such as the quantum 
Hall effect^ and the quantum spin Hall effect, — It is important 
that the system is described by the Dirac Lagrangian, where 
the Green's function behaves as 0(l/|fe|) for k — > oo (k is 
the momentum of the Dirac particle), and thus the ultraviolet 
divergence is likely to occur in contrast to the nonrelativistic 
case with the Green's function oc 0(l/k 2 ). In the study on the 
TI by the diagrammatic expansion, therefore, the divergent di- 



agram should be treated carefully. Since the ultraviolet diver- 
gence is inevitable in calculating physical quantities related 
to the electromagnetic response of the TI, it is a fundamental 
issue how to regularize the divergence. 

Another interesting aspect of TM is the quantization of 
physical observables. In diagrammatic studies, TM appears in 
the bare one-loop Feynman diagram in the self energy of the 
electromagnetic U(l) gauge boson. 7 As shown by the Cole- 
man and Hill, the TM obtained in the 1-loop diagram is ex- 
act, and the correction from scalar, spinor, vector and gauge 
fields is irrelevant^ The coefficient of the TM is quantized, 
and the mechanism can be understood from the viewpoint of 
the Berry's phase. In the band insulator, the correspondence 
between the TM and the Berry's phase is explicitly seen from 
the Hall conductance obtained by the Kubo formula»i2ii£ The 
Berry phase represents the topological number and ensures the 
topological robustness of the TM^i 7 - 

The robustness of the TM against impurities, as it is topo- 
logically expected, is crucially important in obtaining experi- 
mental evidence of the TM. This robustness can be understood 
by the Coleman-Hill theorem in the gapped Dirac systems. 
In diagrammatic calculations, after averaging the positions 
of impurities, the impurities can be viewed as the chargeless 
scalar field and thus the resulting diagrams satisfy the Ward- 
Takahashi identity and the momentum conservation. Hence, 
as long as the pure system free from the impurities is insulat- 
ing and the fermion gap opens, the Coleman-Hill theorem for 
scalar fields can be applied and the TM persists against impu- 
rities. Recently, Nomura and Nagaosa studied the effect of the 
various type of the impurities which can break the time rever- 
sal symmetry and can alter the mass gap on the surface state 
of the TI.— They confirmed that the impurities does not af- 
fect the TM within the linear response theory. Moreover, they 
suggested that the metallic state is localized by the magnetic 
impurities and, even in that case, the magnetoelectric response 
can be described by the TM. 

The gapped surface state of the TI where the TM charac- 
terizes the electromagnetic response is realized by depositing 
the ferromagnet (FM) on the surface of the Tl-si In this sys- 
tem, the magnetization of the FM plays the same role as the 
(7(1) gauge field and the Dirac mass. By using this property, 
various spintronics phenomena have been predicted.—^ For 
example, it has been proposed that the current applied on the 
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surface of the TI can flip the magnetization] 19 - 23 i 27 Nomura 
and Nagaosa showed that the magnetic textures of FM show a 
nontrivial dynamics when the electric field is applied, such as 
the motion of the magnetic domain wall.?- These new effects 
are expected to lead to new spintronics devices. 

In this paper, we demonstrate that the surface of the TI ex- 
hibits the rectification effect because of the parity anomaly 
when the current is pumped by the magnetization dynam- 
ics of FM (see Fig. [TJ. The rectification effect occurs in 
the adiabatic limit, i.e., when the dynamics of the magneti- 
zation is slow, and the validity of the adiabatic approxima- 
tion is studied diagrammatically. To avoid the divergence, we 
adopt the dimensional regularization developed by 't Hooft 
and Veltman.— The dimensional regularization is one of the 
simplest regulators which ensures the gauge invariance of the 
theory, and is generally applicable. In fact, this regularization 
is successfully applied to the study on the conductivity of the 
graphene.— 

We also study the effect of the impurities by using the 
perturbation theory when the energy gap of the surface state 
opens. Although the robustness of the TM is ensured by the 
Coleman-Hill theorem, it may be important to see the robust- 
ness by an explicit diagrammatic calculation. It turns out that 
the first-order result gives the correct quantized TM, regard- 
less of the impurity concentration. 




FIG. 1 : (color online) A ferromagnet (FM) is attached on the surface 
of the TI. The direction of the magnetization M of the FM is exter- 
nally controlled. The current is pumped on the surface of TI by the 
rotation of the magnetization. The precession axis is assumed in the 
y-z plane. 



II. RECTIFICATION EFFECT 

In this section, we study the current driven by the magneti- 
zation precession in a junction of TI and FM. The surface state 
of the TI is described by the (2+l)D massless Dirac Hamilto- 
nian. When a FM is attached, the spin-exchange interaction 



appears and the Hamiltonian is given by 5 

H = J d 2 xw F V/(x)(p y cr K -p x CT y )V'(x)+AV' t (x)(M-cr)V'(x) , 

(2) 

where vp is the Fermi velocity, x = (x, y), cr = (cr x , a y , a z ), 
the direction of the magnetization is represented by a unit vec- 
tor M = (M x , My, M z ) and A is the spin exchange coupling 
energy between the surface fermion and the local magnetiza- 
tion of the ferromagnet. We ignore the effect of the electoro- 
magnetic field^ resulting from the magnetization precession 
since the magnetic field itself does not affect the charge cur- 
rent, and the electric field, which may be proportional to a 
system size^! can be neglected in a small-size system. We 
assume that the Fermi level is at the Dirac point. Then, the 
Lagrangian is simplified as 
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where \x = 0, 1, 2 is the space-time index and 
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= (t,x), = diag(vf\-l, -1) 



Here, 7 M are the 2x2 matrices and satisfy the Dirac alge- 
bra {7^, = 2g'" y x 1-2x2 (1-2x2 is the 2 x 2 unit matrix). 
The Lagrangian is thus equivalent to the conventional Dirac 
fermion of the charge e(< 0) coupled with the electromag- 
netic field. In the early 1980's, it was found that the current is 
described by the TMi£: 
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where = vpeip(x)^f ,1 ij;(x), the summation over v, p is im- 
plied and |ft<9iA/A|, \hv^ l dtA/A\ <§; mvp. In terms of the 



magnetization, the current is given as discussed in Ref . 1251 by 



(j°(x))=p(x) 
(f (x)) ■■ 
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where i = (x, y). 

These results are for the case of constant to, but they re- 
main correct as long as the adiabatic condition is satisfied, 
namely, if the mass of the Dirac particle, to, is larger than the 
frequency of the magnetization multiplied by Hluq The va- 
lidity of eq. © in the adiabatic regime will be demonstrated 
in Sec. [Ill] and we will proceed here based on eq. ©. From 
eq. ©, we see that and (f ) exhibit anomalous behav- 
ior when to crosses zero because of the coefficient to/|to| = 
sign(TO). In fact, if dtM is finite when to = 0, (J 1 ) exhibits 
a sudden jump. If dtM changes sign at m = 0, the sign of 
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the current is the same for both m > and m < 0, and the 
generated current is rectified due to the anomaly. 

Let us look into these behaviors in detail in the case of pre- 
cessing magnetization with the angular frequency ljq. We 
choose the angle (f> as the angle of the precession axis mea- 
sured from the z-axis (Fig.|TJ, and the radius of the precession 
is denoted by ro(< 1). The magnetization is then represented 
as 



r sin LUot 
M = | ro cos cf> cos w i + y/T- 
- ro sin <j> cos wot + yl 



(7) 



Tq cos < 



The induced current is given by — and 
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-ro cos 4> sin woi . 



Let us examine the case of 
magnetization reads 



7r/2 in detail, where the 
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(ro sinw t, \Jl-r%, -r cos ut) 



(9) 



The mass term is thus m = — roA cos wot/vp and eq. © re- 
duces to = and 



(f) 



eA 
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(10) 



In this case, the generated current is a direct current as shown 
by the straight line in Fig. [2] in spite of the oscillation of the 
sign of the magnetization. By using the reasonable material 
parameters 2 ^ as — 4 x 10 5 m/s, A = 0.2 eVr& the cur- 
rent density is given by (j x ) = 9.6 x 10 _12 rowo| cos u>ot\ 
A/m. If the frequency w is 1GHz, (j x ) = ~ 
A/m, which is large enough for detection. When (f> devi- 
ates from it/2, i.e., <f> = it/2 + 6<fr, (j x ) exhibit a jump 
Jtanty (= Tf^v/l^cot^) 



between ±|^^ 



when m = 0, namely, when rocoswot = — yl — rg tani50. 
A jump becomes significant when ro ~ 1/V2 and <fi is close 
to 7r/4 as shown by the dashed line in Fig. [3] 



III. VALIDITY OF ADIABATIC APPROXIMATION 

We have shown above that a topological effect of the parity 
anomaly is directly observable either as a rectification effect 
or as a jump in the electric current. As we have pointed out, 
however, eq. © is not correct in the vicinity of m = due 
to the non-adiabaticity. In this section, we investigate the ef- 
fect of non-adiabaticity. The result of eq. © is correct as far 
as wo -C mvp / h (adiabatic limit) since the gradient expan- 
sion is implied. In the opposite case of mVp/H < ujq, the 
time-dependence of the "mass" term becomes essential and 
a deviation from the adiabatic limit (eq. ©) may arise. We 
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FIG. 2: (color online) The current on the surface of the TI induced by 
the rotating exchange field for <f> = it/2, tv/2 + 0.2 and ro = l/\/2- 
The straight line denotes the full-wave rectification given by eq. (TO). 
The dashed line represents the current when the rotation axis of the 
exchange field deviates from x-y plane. The abrupt change occurs at 
M z = and will be smoothed by the high frequency mode in the 
short time period St ~ h/ A. 
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FIG. 3: (color online) The current for <f> = 7r/4, 7r/4 + 0.05, and 
ro = 1 /V2. For cf> — n/4, the z-component of the rotating exchange 
field M z touches zero. For <f> > n/4 the abrupt change of the current 
appears when M~ = 0. 
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will demonstrate here explicitly that the non-adiabatic correc- 
tion from the dynamic mass term (m(t)) is negligibly small 
and that the time window for the non-adiabatic correction 
to be dominant is very narrow in the low frequency regime 
uiq <C A/H, which practically justifies the adiabatic treatment. 
For A = 0.2 eV, this condition is validated for ljq <C 400 
THz. 

In the following calculation, we set up = 1 and H = 1 for 
simplicity. To take into account the effect of the dynamical 
mass, we introduce the small fluctuation, m 2 (a;), as 



and 



m{x) = to + 7712(2;) 



(11) 



where m = 7^7 J d 3 xM z {x) is the average mass (TV is the 
volume of the spacetime) and x = (t, x) denotes space-time 
coordinates. The mass fluctuation leads to an additional po- 
tential term in the Lagrangian: 



SL = —m2(x)ip(x)4>(x) ■ 



(12) 



Including the correction to the linear order in m 2 (Fig. the 
current becomes 



W(k))=iA v (-k) 



where 



-im 2 (-fc) 



S{k) 



d 3 p 
(2tt) 3 
d 3 p 



tt[S(k + P )S(p)-f] , 
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k 2 — to 2 

and k 2 = (k ) 2 - (fc 1 ) 2 - (fc 2 ) 2 . It is calculated as 



n^(fc) = 



mfc^ 1 



dx- 



(1 - 2x) 



it j (m 2 
( = for k 2 < 4m 2 ) 



x(i - x)k 2 y/ 2 



(17) 



The current within the linear response theory is obtained by 
substituting fc° — > k° + i0 + . In this calculation, the ultraviolet 
divergence is regularized^ - by the dimensional regularization 
of 't Hooft and Veltman^ The dimensional regularization is 
useful in the sense that the obtained result automatically satis- 
fies the gauge invariance of the theory. To regularize the diver- 
gent diagram, there are various types of regulators other than 
the dimensional regularization. For example, the hardcutoff 
that assumes the momentum cutoff in the divergent integral is 
a simple regulator. However, as known in the gauge theory 
of the Dirac fermion, the resulting current from the hardcutoff 
does not satisfy the Ward-Takahashi identity. — 

As seen from eqs. ( fl~6b . (T7\ . TM term is dominant when 



fc" <C m . In the time-dependent system we consider, 
fc 2 ~ (fc ) 2 = lj 2 , and thus the adiabatic condition where 
eq. (|6]i is correct is given by uj <C m. When m ~ 0(u!q), 

(13) all the contributions in eq. (TT~5T > becomes the same order, 
111 ~ ri2 P ~ II ~ 0((jJq) as we see from dimensional anal- 
ysis. However, the time-window for m < ojq is H/A and 
is very narrow if luq <C A/H. The non-adiabatic contribu- 
tion is therefore essentially not observable if ljq <C A/h, and 

(14) then the rectification effect and sudden jump in the pumped 
current due to the anomaly would be observable even in the 
non-adiabatic regime. 





FIG. 4: The diagrams that contribute the current in the first order in 
A and ni2. 



IV. EFFECT OF IMPURITIES 

Next, let us consider the effect of impurities on the massive 
Dirac fermion. The TM is topologically protected and robust 
against impurities in the gapped phase in contrast to the metal- 
lic phase where impurities induce a nontrivial dynamics. 32 
The effect of the impurities in the Hamiltonian is given by 



(f{k)) = or 



where 39 



fc^fc" 



)n 1 {k) + e^ v n 2p {k))A 1/ (-k) 



+ H^(fc)m 2 (-fc), 



(15) 
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(16b) 



Hu 



d 2 x?j i (x)^ t (x)V'(x) , 



(18) 



where Vj is the impurity potential. In this paper, we con- 
sider fi(x) = U( ^( x — x n)' where u is the strength of 
the impurity scatterers, Ni is the number of the impurities 
and x n is the position of the n-th impurity. We will aver- 
age over the positions of the impurities as (vi(x))i = and 
(i>i(xi)i!i(x 2 ))i ~ S(xi — x 2 ), where V is the volume 
of the space. For treating the dissipation driven by the impu- 
rities, we use the Keldysh Green's function. The current is 
given by 



(j l (x)} = iehvpEijtila^ G < (x , x)} , 



(19) 



where i, j = (x,y) and G < is the lesser Green's function, 
G< a ,(x,x') = I^'OeOVvOe)). In the first order in A, the 
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current reads 

x[/(fco + y)-/(fco-y)] 

x tt[o*g r (k - |) A a (fc + §)] , (20) 

where / is the Fermi distribution function. Here the advanced 
(retarded) Green's function, g a [= (g r )*], is given as 

_ (fefcp - z??)l 2 x2 + m^tr* - hv F (k x a y - k y a x ) 
9 W ~ {hk a - irj) 2 - (mvlY - (^ F |k|) 2 

(21) 

where r} represents the self energy of the fermion.— To see 
the effect of the impurities on TM, we focus on the term pro- 
portional to tr[er'' <j v ' a p \ = 2ie^ vp . As a result, it is found 
that TM term does not change even if the impurities exist. 
Now, let us discuss the effect of the higher order diagrams 
in the context of the Coleman-Hill theorem. 11 This theorem 
states that the correction from scalar, spinor, vector and gauge 
fields does not affect TM in the photon self energy appearing 
in the 1-loop diagram, if the Ward-Takahashi identity is sat- 
isfied, the momentum is conserved and the Green's function 
is massive. By averaging over the positions of impurities, the 
impurity can be viewed as the chargeless scalar field, and then 
the Ward-Takahashi identity and the momentum conservation 
is satisfied. Because of the averaging process, the impurity by 
itself does not induce the current ((i>j(x)), = 0), in contrast 
to the 777,2 term discussed in the previous section. Therefore, 
the Coleman-Hill theorem 11 can be applied to our calculation 
and the impurities do not affect TM at any order of loops. 

V. CONCLUSION AND DISCUSSION 

We have discussed the current dynamics on the surface of 
the TI by the precessing magnetization. The rectification ef- 
fect and the jump of the current have been found to arise from 
the parity anomaly in the adiabatic regime. We have discussed 
that the effect of the non-adiabaticity can be neglected. To 
remove the divergence in the diagram preserving the gauge 
invariance of the theory, the dimensional regularization2£ is 
used. It should be pointed out that the dimensional regular- 
ization can be also applied to the calculation of other physical 
quantities such as the conductivity in the metallic phase of the 
TI and hence would be a useful tool to investigate properties 



of TI. Also, the robustness of the TM against the impurities is 
reconfirmed. 

Finally, let us discuss the possibility of a rectification driven 
by a light. In our study, it is essential for the rectification 
that the mass term controls the sign of the Hall current. The 
same goes for the case of a circularly polarized light radia- 
tion to the surface of a bare TI, and thus a rectification is ex- 
pected. In this case, the light produces a varying mass term 
777 = g/j,BB z (x)/vp by the Zeeman coupling, where B z is 
the magnetic field of light. The Zeeman couplings to B x , B y 
does not produce an observable current since these are much 
weaker than the exchange coupling discussed in Sec. [TT] On 
the other hand, the rotating electric field of light induces the 

2 

current = §j;-ffieijEj with i,j = (x,y). If the pointing 
vector is within the surface plane (B x E) z = 0, E x and E y 
change the sign when B z = 0; the rectification occurs. For 
experimental realization, a strong laser would be necessary to 
satisfy the adiabatic condition Hui, hvpk <C g^s\B z \. In addi- 
tion, the fine tuning of the Fermi level at the Dirac point needs 
to be carried out precisely, since the induced gap is assumed 
to be small, for example, 10 -5 eV for B z = 100 mT. 

We expect that the nontrivial-current dynamics driven by 
the quantum anomaly can be applied to new spintronics de- 
vices. 

Note added. After the completion of our work, we became 
aware of the paper by B. Dora et ai, who studied a rectifi- 
cation effect on a edge state of 2D TI induced by a circu- 
larly polarized light both in the adiabatic and non-adiabatic 
regime.-- The edge state is described by the (1+1)D-Dirac 
fermion. In their setting, the rectification in the adiabatic 
regime is driven by the topological excitation 5 ^ related to the 
Dirac mass term. On the other hand, in our (2+l)D study, the 
parity anomaly plays an essential role in the rectification. 
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